To APPEAR. IN THE JANUARY 2002 ISSUE OF ApJ (VOL. 565) 

Preprint typeset using lAT^^X style cmulatcapj 



EXPLOSIVE NUCLEOSYNTHESIS IN ASPHERICAL HYPERNOVA EXPLOSIONS 
AND LATE TIME SPECTRA OF SN1998BW 

Keiichi Maeda^, Takayoshi Nakamura^, Ken'ichi Nomoto^'^, Paolo A. Mazzali^'^, Ferdinando Patat^, Izumi 

Hachisu"' 

To appear in the January 2002 issue of ApJ (Vol. 565) 



o : 
o . 

(N ■ 

^• 
(D ■ 

c/^ ; 
i> ■ 

(N : 

(N ■ 

> • 

o 

o ■ 



o . 
o ■ 

Oh 
I 

O ■ 



ABSTRACT 

Aspherical explosion models for the hypernova (hyper-energetic supernova) SN 1998bw are presented. 
Nucleosynthesis in aspherical explosions is examined with a two-dimensional hydrodynamical code and a 
detailed nuclear reaction network. Aspherical explosions lead to a strong a-rich freezeout, thus enhancing 
the abundance ratios ['*^Ca, ''^Ti, and ^^Zn / Fe] in the ejecta. The nebular line profiles of the Fe- 
dominated blend near 5200 A and of [O i] 6300,6363 A are calculated and compared with the observed 
late time spectra of SN 1998bw. Compared with the spherical model, the unusual features of the observed 
nebular spectra can be better explained if SN 1998bw is a strongly aspherical explosion with a kinetic 
energy of ~ 10^^ ergs viewed from near the jet direction. 

Subject headings: supernovae: individual: SN1998bw — gamma rays: bursts — nucleosynthesis — line: 
profiles 



1. INTRODUCTION 

The exceptionally bright Type Ic supernova (SN Ic) SN 
1998bw was discovered as the probable optical counterpart 
of the gamma-ray burst GRB980425 (Galama et al. 1998). 
The early light curve and the spectra of SN 1998bw have 
been successfully modeled as the hyper-energetic explosion 
(kinetic energy i? ~ 4 x lO'^^ ergs) of a massive C-t-0 star 
(Iwamoto et al. 1998; Woosley et al. 1999; Branch 2001). 
In this paper the term "hypernova" is used to refer to a 
SN explosion with E > 10^^ ergs, regardless of the nature 
of the central engine (Nomoto et al. 2001a). 

Despite the success of the hypernova model in repro- 
ducing the observed features of SN 1998bw at early times, 
some properties of the observed light curve and spectra at 
late times are difficult to explain. (1) The tail of the ob- 
served light curve declines more slowly than the synthetic 
curve, indicating that at advanced epochs 7-ray trapping 
is more efficient than expected (Nakamura et al. 2001a; 
SoUerman et al. 2000). (2) In the nebular epoch, the 
[O i] 6300 A emission is narrower than the emission near 
52OOA. As discussed in Mazzah et al. (2001), this latter 
feature is mostly due to a blend of [Fe 11] lines. Mazzali 
et al. (2001) calculated synthetic nebular-phase spectra of 
SN 1998bw using a spherically symmetric NLTE nebular 
code based on the deposition of 7-rays from ^^Co decay 
in a nebula of uniform density and composition. They 
showed that the [O i] and the [Fe 11] features can only 
be reproduced if different velocities are assumed for the 
two elements. A significant amount of slowly-moving O is 
therefore necessary to explain the profile of the [O i] line. 

Both these features are in conflict with what is expected 
from a spherically symmetric explosion model, where 7-ray 



deposition decreases with time and where iron is produced 
in the deepest layers and thus has a lower average velocity 
than oxygen. Mazzali et al. (2001) suggested that these 
are signatures of asymmetry in the ejecta. Therefore in 
this paper we examine aspherical explosion models for hy- 
pernovae. 

Aspherical explosions of massive stars have been inves- 
tigated as possible sources of gamma-ray bursts (GRBs) 
(Woosley 1993; Paczynski 1998). MacFadyen & Woosley 

(1999) showed numerically that the collapse of a rotat- 
ing massive core can form a black hole with an accretion 
disk, while a jet emerges along the rotation axis. The jet 
produces a highly asymmetric explosion (Khokhlov et al. 
1999). However, these studies did not calculate explo- 
sive nucleosynthesis, nor did they show spectroscopic and 
photometric features of aspherical explosions. Nagataki 

(2000) performed nucleosynthesis calculations for aspher- 
ical SN explosions to explain some features of SN 1987A, 
but he only addressed the case of a normal explosion en- 
ergy. 

In the present study, we examine the effect of aspheri- 
cal explosions on nucleosynthesis in hypernovae. We then 
investigate the degree of asphericity in the ejecta of SN 
1998bw, which is critically important information to con- 
firm the SN/GRB connection, by computing synthetic 
spectra for the various models viewed with different orien- 
tations and comparing the results with the observed late 
time spectra of SN 1998bw (Patat et al. 2001). 

2. ASYMMETRIC EXPLOSION MODELS 



The first step of our calculation is the hydrodynamical 
simulation of the explosion with a 2D Eulerian hydrody- 
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namical code based on Roe's scheme (Hachisu et al. 1992, 
1994). Eulcr's equations arc solved with a constant adi- 
abatic index 7 = 4/3, which is a good approximation if 
the pressure is radiation-dominated. The effect of nuclear 
reactions on the hydrodynamics is negligible since the ex- 
plosion energy is large. We use 120 x 120 meshes on a cylin- 
drical (r, z) coordinate system. The; nicsh size is linearly 
zoned and decreases inward, which gives a high resolution 
of the hydrodynamic evolution of the central regions where 
explosive nucleosynthesis takes place. We follow 190 test 
particles initially in the Si layer and 2250 particles initially 
in the C-l-0 layer, tracking their density and temperature 
histories. These histories are then used to calculate the 
change in the chemical composition, using a reaction net- 
work including 222 isotopes up to ^^Ge (Thielemann et al. 
1996). 

We construct several asymmetric explosion models for 
various combinations of the model parameters (Table 3). 
We use as progenitor the 16 Mq He core of a 40 Mq star 
(Nomoto & Hashimoto 1988). This has a 13.8 M© C+0 
core, the same as that used in Iwamoto et al. (1998). We 
test three values of the final kinetic energy: E = lx 10^^ , 
1 X 10^^, and 3 x 10^^ ergs. The hydrodynamical simu- 
lation is started by depositing the energy below the mass 
cut that divides the ejecta from the collapsing core. The 
energy deposited is divided between thermal and kinetic 
energy, with various ratios. The asymmetry is generated 
by distributing the initial kinetic energy in an axisymmct- 
ric way. This is done by imposing different initial velocities 
in difFcrcnt directions: Vz = a z in the jet direction and 
Vr = P r on the equatorial plane. The ratio a/ (3 ranges 
from 16:1 to 1:1 (spherical case). The mass cut is set at 
Mr = 2AMq, so that the ejected mass of ^'^Ni is ^ 0.4 
Mq to reproduce the peak of the light curve (Nakamura 
et al. 2001a) in models A, B, C, E, and F. 

3. NUCLEOSYNTHESIS 

Figure 1 and 2 show respectively the post-shock peak 
temperatures and densities for the asymmetric hypernova 
model C in the direction of the jet and perpendicular to it 
(with those for the spherically symmetric hypernova model 
F also shown for comparison), and the isotopic composi- 
tion of the ejecta of model C. In the z-direction, where 
the ejecta carry more kinetic energy, the shock is stronger 
and post-shock temperatures are higher, so that explo- 
sive nucleosynthesis takes place in more extended, lower 
density regions compared with the r-direction. Therefore, 
larger amounts of a-rich freeze-out elements, such as "^He 
and ^^Ni (which decays into ^^Fe via ^^Co) are produced 
in the z-direction than in the r-direction. Also, the ex- 
pansion velocity of newly synthesized heavy elements is 
much higher in the z-direction. The velocity of elements 
ejected in the z-direction in model C is actually similar to 
the result of a spherical explosion with 3 x 10^^ ergs 
(Nakamura et al. 2001b), although the integrated kinetic 
energy is only E = 1 x 10^^ ergs. 

In contrast, along the r-direction ^^Ni is produced only 
in the deepest layers, and elements ejected in this direc- 
tion are mostly the product of hydrostatic nuclear burning 
(O), with some explosive oxygen-burning products (Si, S, 
etc). The expansion velocities are much lower than in the 
^;-direction. 



Figure 3 shows the 2D distribution of Ni and O in 
model C in the homologous expansion phase. Near the 
z-axis the shock is stronger and a low density, "'He-rich re- 
gion is produced. ^^Ni is distributed preferentially in this 
direction, but it is mostly located slightly off of it because 
the shock propagates laterally as it penetrates the stellar 
envelope. As a result, the distribution of heavy elements 
is elongated in the z-direction, while that of ^^O is less 
aspherical. On the other hand, because the ejecta move 
more slowly in the r-direction, densities in this direction 
are higher than in the z-direction. 

Tables 1 and 2 give respectively the detailed yields and 
the abundances of major stable isotopes relative to the so- 
lar values for model C (in Table 2, [A/B] = log^o (A/B) - 
logj^Q (A/B)q, where A and B are nuclear mass fractions). 
The main characteristics can be summarized as follows (see 
also Nomoto et al. 2001b). 

(1) The complete Si-burning region is more extended for 
larger explosion energies. The aspherical explosion causes 
a region of higher entropy along the z— axis, which offers 
better conditions for the a-rich freezeout (Fig. 1). The 
high entropy inhibits the production of ^^Ni. Much ^He 
is left after the freezeout, so that the elements produced 
through ^He capture are very abundant in the deepest 
region along the z-axis (Fig. 2). This residts in the en- 
hancement of the elements synthesized in the deepest re- 
gion, such as *^Ca (produced as ^*Ti), ''^Ti (as ^*Cr), and 
elements heavier than A ^ 58. Because of the enhance- 
ment of these elements and the simultaneous suppression 
of ^^Ni, the abmidances of these elements relative to iron 
(e.g., [^''Ca, ^^Ti, ^'^Zn /Fe]) are greatly enhanced. For 
more asymmetric explosion, the effect of a-rich freezeoiit 
is even larger. 

(2) Incomplete Si-burning and 0-burning regions are 
more extended for larger explosion energies (Nakamura et 
al. 2001b). This resuhs in the enhancement of ^^Si, ^^S, 
^°Ca, ^^Cr (produced as ^^Fe), ^''Fe, and in the reduc- 
tion of O. Asphericity has little effect on the production 
of these elements. 

The most pronounced effect of asphericity is that ele- 
ments produced by the strong a-rich freezeout are greatly 
enhanced relative to iron (e.g., [Ti/Fe]). For other ex- 
plosive burning products, the effect of a large explosion 
energy usually dominates over that of asphericity. 



4. THE LATE TIME SPECTRA OF SN 1998BW 

In order to verify the observable consequences of an 
axisymmetric explosion, we calculated the profiles of the 
[Fe 11] blend and of [O i] for models A-G. Line emissivi- 
ties were obtained from a ID NLTE nebular code (Mazzali 
et al. 2001), and the column densities of the various ele- 
ments along different lines of sight were derived from the 
element distribution obtained from our 2D explosion mod- 
els. Because we assume that the nebula is optically thin, 
the blended nature of the emissions is automatically taken 
into account. 

These are compared to the 26 Nov 1998 spectrum of 
SN 1998bw. We select this spectrum because the wave- 
length of the 52OOA feature, which was somewhat redder 
at earlier epochs, at this and later epochs coincides with 
that of the equivalent [Fe 11] feature in the nebular spectra 
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of SNc la (Axclrod 1980), indicating that other contri- 
butions (sec Section 5) are now neghgiblc. Tabic 3 gives 
the FWHM of our synthetic hnes as a function of viewing 
angle. The corresponding observed value is 380A for the 
Fc-blcnd. This is estimated assuming that the continuum 
level is the value around 5700A where the flux has the 
minimum value. 

The FWHM of the 5200A feature is narrower at this 
epoch than at earlier ones, and so is that of [O i]. While 
for the former feature other contributions may be respon- 
sible for a broader line at earlier epochs, in the case of [O i] 
the decreasing density of the outer envelope must be the 
principal reason. At late epochs the density of the outer 
envelope is expected to be too small to trap the 7-rays. 
Therefore wc set the oiitcrmost velocity of the emitting 
region to reproduce the FWHM of the [O i] line (150A), 
and then calculate the profiles of the [Fe 11] blend. 

The [Fe 11] and [O i] profiles for model C viewed at an 
angle of 15° from the jet direction and those for model F 
are compared to the observed spectrum of 26 Nov 1998 in 
Figure 4. For model F in Fig. 4, the outermost velocity of 
the emitting region is set to make the [Fe 11] line as broad 
as possible, because for this model wc cannot get a rea- 
sonable fit for the [O i] line, which is always much broader 
than the observation. Indeed, fitting the [O l] line was not 
possible for all models. Among the hypernova models, in a 
spherical explosion (model F) oxygen is located at higher 
velocities than iron, and the [O l] line is too broad for any 
choice of the outer velocity of the emitting region. This 
is because of the deficiency of oxygen with small velocity 
along the line of sight. Also, even though the Fe feature 
can be wider than the O line if O and Fe are mixed exten- 
sively, the expected ratio of the width of the Fe-blend and 
the O line even in a fully mixed model is ~ 3 : 2 (Mazzali 
et al. 2001). This is the result of taking blending into ac- 
count, but giving all contributing lines the same intrinsic 
width. However, the observed ratio is even larger, > 2 : 1, 
implying that the [Fe 11] lines are intrinsic; ally broader than 
the [O i] ones. Therefore the observed line profiles are not 
explained with a spherical hypernova model. The same is 
true for the moderately asymmetric model E viewed near 
the equator. 

In our aspherical explosion models Fe is distributed pref- 
erentially along the jet direction, and so a larger ratio of 
the Fe and O line widths can be obtained. All the strongly 
aspherical hypernova models A, B and C, when viewed 
from a near-jet direction, give line widths comparable to 
the observed values. The very energetic model B cannot 
reproduce the O line when viewed near the equator, but 
this is because O is too fast near the equator and too de- 
pleted near the poles to give a low-velocity component. 

When the degree of asphericity is high and the explosion 
is viewed from near the jet direction, the component lines 
in the [Fe 11] blend have double-peaked profiles, the blue- 
and red-shifted peaks corresponding to matter situated in 
the two opposite lobes of the jet, where Fe is mostly pro- 
duced. Because of the high velocity of Fe, the peaks are 
widely separated, and the blend is wide (Fig. 4, model C). 
In contrast, the [O l] line is narrower and has a sharper 
peak, because O is produced mostly in the r-direction, at 
lower velocities and with a less aspherical distribution. 

Figure 5 shows the [O i] line for model C at different 



orientations. The mean expansion velocity of O is lower in 
the aspherical cases than in the spherical model, because in 
aspherical models low-velocity, high-density 0-dominated 
matter is found near the center (Fig. 3). Therefore the 
width of the O line in the aspherical models (A, C) viewed 
from near the equator can also be comparable to the ob- 
served width. This, however, does not mean that the line 
profile is always sharp as seen in SN1998bw. As shown 
in Figure 5, when the angle is large, the line first broad- 
ens, and eventually it develops two peaks. The reason 
can be seen in Figure 3. The highest density region has 
a typical velocity ~ 3000 km s~^ along the r-axis in this 
model. This corresponds to a Doppler shift of ^ 120A 
between the approaching and receding parts when the SN 
is viewed from the r-direction. Therefore, to produce the 
narrow and sharply peaked O line in a hypernova model, 
the explosion must be aspherical and viewed from near the 
polar direction. 

Figure 6 shows the profiles of the [Fe 11] blend viewed 
at 5° for various aspherical hypernova models (A, B, C 
and E). In all of these models, the computed O line re- 
produces the observed one as discussed above. The Fe- 
blends appear to reproduce the 5200A feature reasonably 
well. The profiles shows small peaks, which arc not seen 
in the observations. These peaks are more pronounced 
in more asymmetric models. Mixing of the ejecta may 
distribute the ^^Ni to lower velocities, thus reducing the 
double-peaked profiles of the Fe lines. A spherical hyper- 
nova model (model F in Fig. 4) also gives a broad Fe line, 
and without the sharp peaks. However the O line is much 
too broad. Also, in this work we have used a spherical, 
uniform density nebular model to compute line emissiv- 
ities. In an aspherical model, dense central and equa- 
torial regions may have higher 7-ray trapping efficiency, 
which may result in stronger low-velocity line emission 
than in our model. Thus the component iron lines in the 
blend may have wider, flat-top proflles rather than double- 
peaked shapes, which coiild eliminate the minor peaks in 
the Fe blend seen in our present models. 2D 7-ray trapping 
calculations are therefore needed to compute the detailed 
spectra and the light curve. 

It is reasonable to think that asphericity reduces the 
energy below that estimated previously based on spheri- 
cally symmetric models. To examine this, we turn now to 
the lower energy explosion models (D, G). The [Fe 11] and 
[O i] profiles for these models are shown in Figure 7. First 
of all, these models always give a narrower O line than 
the observed one. In these low energy models, in fact, the 
velocity of the ejecta is too small. The fastest-moving mat- 
ter approaches the observer with a velocity of 3500 km s~^ 
for model G, and 5000 km s^^ in the case of the aspheri- 
cal model (D) viewed from the z-direction. The observed 
FWHM (150A), however, indicates that there should be 
material moving faster than 7000 km s"'^. Also, because 
of the low velocity, the component [Fe 11] lines are too 
narrow and do not blend to form a broad feature. Prom 
these arguments, we conclude that the explosion energy of 
SN1998bw should have been large, Ek ^ 1 x lO^^ergs. 

5. CONCLUSION AND DISCUSSION 

We calculated the nucleosynthesis in aspherical hyper- 
nova explosions. We find that in such explosions Fe is 
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mostly ejected at high velocity in a jet along the polar di- 
rection, while nearer the equatorial plane burning is less 
effective, and low-velocity O is mostly ejected. We show 
that the unusual ratio of the width of the O and Fe neb- 
ular lines in SN 1998bw can be explained with a strongly 
aspherical explosion model viewed from a near-jet direc- 
tion. Also, in this case the O line has a very sharp peak, 
in agreement with the observations. 

Much of our discussion was based on the identification 
of the 52OOA feature as a blend of [Fe 11] lines. Although 
several caveats apply to that identification, because other 
lines may contribute, we claim that [Fe 11] lines dominate. 
Other possible contributions are as follows. According to 
Sollerman et al. (2000) the feature contains not only [Fe 11] 
lines but also lines of [Mg i] , [O i] and possibly [Fe i] . Lines 
of [Mg l] are included in our nebular code, and the strength 
of the 5 172 A emission, which is consistent with the 5470A 
line reproducing the observed [Mg i] peak near 45OOA IS 
much smaller than that of the [Fe ii] lines. [O i] 5577A is 
also included in our code, and it is strong and fills up the 
red part of the emission, but it does not contribute to the 
blue side (Mazzali et al. 2001). As for [Fe i], Sollerman 
et al. (2000) say that their models 'have too low a degree 
of ionization', suggest that the density in the Fe-emitting 
regions should be lower, and 'regard the Fel emission as 
dubious'. We confirm this result. 

Another possibility is allowed Fe II emission, such as 
identified by Fihppenko (1989) in the Type II SNe 1987F 
and 19881 and by Filippenko et al. (1990) in the Type 
Ic SNe 1985F and 1987M. The main Fe 11 features are at 
4570A (muhiplets 37 and 38), 5190A (multiplet 42) and 
5320 A (multiplets 42 and 43). These lines are not in- 
cluded in our atomic model, because not all the collision 
strengths for these transitions are available. In particular 
multiplet 42 is not available. However, several arguments 
apply against the Fe 11 identification in SN 1998bw. 1) In 
SNe 1987F and 19881 all three features are strong, while in 
the SNe 1985F and 1987M the feature at 4570A is narrow 
and it is identified as [Mg l], suggesting that the feature 
at 52OOA could be [Fe 11]. A similar situation holds for 
SN 1998bw, also a SN Ic. 2) Fe 11 emission occurs at 
high density. Filippenko (1989) infers log rig ~ 9 - 10 for 
SN 1987F. However, the densities we derive for SN 1998bw 
are more than one order of magnitude smaller (Mazzali 
et al. 2001). 3) The relatively low density is confirmed 
by the large ratio (~ 4) of the Ca 11 IR emission com- 
pared to [Ca 11] 73OOA. Fig.2 of Ferland & Persson (1989) 
suggests that log rig < 9- Finally, the models of Mazzali 
et al. (2001), which are based on [Fe 11] only, reproduce 
the 52OOA feature very well. In any case, extending the 
model to include Fe 11 lines is a worthwhile effort which 
we arc going to make. 

The 52OOA feature is bluer on Nov. 26 and thereafter 
than on earlier epochs. Two principal factors are probably 
responsible for that: 1) the disappearance of the contin- 
uum, which is still significant in the earlier spectra; 2) the 
reduced intensity of O i] 5577A, as shown also in our mod- 
els (Mazzali et al. 2001). This is a high-density line, and 
it decreases quickly in strength as the ejecta expand. Fad- 
ing Fe II emission may also play a role in causing a small 
wavelength shift. 

Weak unaccounted emission is present to the blue of the 



[Fe 11] feature (4800-5100A). A similar emission is present 
in the nebular spectra of SNe la as well (Patat et al. 2001), 
and it is poorly reproduced by synthetic spectra. It may 
be due to a forest of weak [Fe 11] and [Co 11] transitions 
whose atomic parameters arc not well known. 

The [O i] line in other SNe Ib/c (Filippenko et al. 1995, 
Matheson et al. 2001) also shows a strongly peaked pro- 
file, as in SN 1998bw. This probably signals the existence 
of oxygen at low velocity in most SNe Ib/c. Although 
there are very few SNe Ib/c in which the feature centered 
at 5200 A is detected with enough S/N, it appears that in 
these objects the [O i] line has a similar profile. For low 
energy supernovae (1 x 10'''^ergs), the spherical model G 
shows a sharply peaked [O i] line profile (Fig. 7). This 
seems to favor spherical explosions for these low enc;rgy 
supernovae. However, it is premature to conclude that 
asymmetry is completely absent in low energy SNIb/c, be- 
cause the effects of asphericity in a late time spectra is not 
so large in case of normal supernovae as seen in Figure 7. 
This is due to the low expansion velocities. The oxygen 
line can be narrow even in the case of an aspherical ex- 
plosion viewed from near the equator. For example, it is 
difficult to distinguish model D viewed from 75 deg from 
model G. Because the probability of viewing an asphericl 
SN from the near polar direction is smaller than it is to 
view it near the equator, the observations are not incon- 
sistent with the possibility that most SNe Ib/c are more 
or less aspherical. More nebular spectra of SNe Ib/c with 
higher S/N ratio are needed for further investigation. 

Our aspherical explosion models may be able to explain 
the slow decline of the late light curve of SN 1998bw. In 
these models, the; equatorial region is denser than in a 
spherically symmetric model with the same explosion en- 
ergy. At advanced epochs this region may be able to trap 
7-rays more efficiently than a spherical model, as first sug- 
gested by Nakamura et al. (2001a). Chugai (2000) showed 
that a spherical model could reproduce the light curve, 
if the density near the center (i.e. in the Fe-dominated 
region) was increased above that of the hydro dynamical 
model of Iwamoto et al. (1998). However, Sollerman et al 
(2000) find that the 0-dominated region should be dense, 
and the Fe-dominated region not dense. Although taken 
individually these conclusions are probably correct, they 
appear to be in confiict with one another. An aspher- 
ical model offers a natural solution, because it predicts 
the presence of high density, 0-dominated matter near the 
center (Fig. 3). 

A small degree of linear optical polarization (~ 0.5%) 
was reported in SN 1998bw (Kay et al. 1998; Iwamoto 
et al. 1998; Patat et al. 2001). This can be explained with 
different combinations of asphericity and viewing angle, or 
with large scale clumping in a basically spherical envelope. 
One solution is that ejecta with a moderate departure from 
sphericity are viewed from slightly off the axis of symme- 
try. Our strongly aspherical explosion model has an axis 
ratio of about 3:2 at the outer edge of the oxygen enve- 
lope (Fig. 3), therefore it is consistent with the observed 
polarization, if it is viewed from near the jet direction. 

This work has been supported in part by the grant-in- 
Aid for COE Scientific Research (07CE2002, 12640233) of 
the Ministry of Education, Science, Culture, and Sports 
in Japan. 
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Table 1 

Detailed Yields of Model C (M©) 



Species 


[X/OJ- 


[X/Fe 


Species 


[X/OJ 


[X/FeJ 


1^0 


-1.31 


-0.880 


40Ca 


-0.183 


0.246 




-4.09 


-3.66 


44 Ca 


0.0989 


0.529 


160 


0.000 


0.430 


45Sc 


-1.17 


-0.740 


19p 


-5.40 


-4.97 


48Ti 


0.0153 


0.445 


20Ne 


-0.437 


-0.00731 


51v 


-0.538 


-0.109 


23Na 


-0.210 


0.220 


52 Cr 


-0.215 


0.215 


24Mg 


-0.146 


0.284 


55Mn 


-0.605 


-0.176 


27 Al 


0.127 


0.557 


54Fe 


-0.224 


0.206 


28Si 


-0.0516 


0.378 


56 pg 


-0.430 


0.000 


31p 


-0.0192 


0.411 


59C0 


-0.668 


-0.238 


32s 


-0.182 


0.248 


58Ni 


-0.220 


0.210 


35 CI 


-0.769 


-0.339 


03 CU 


-1.46 


-1.03 


30 Ar 


-0.244 


0.186 


04Zn 


-0.750 


-0.320 


3''K 


-0.979 


-0.549 









Table 2 

Abudances of major stable isotopes relative to the solar values for Model C 



«[A/B] = logio{A/B) - logio(A/B)0 
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Model a 113 




i^K fraction^ 


anele'^fdegl 


Fe FWHM'=[A] 


O FWHM [A] 


outer velocity [km s~^] 


A 16 


10 


0.65 


5 


376 


150 


9000 








15 


344 


150 


8500 








30 


304 


150 


8000 








45 


248 


150 


6500 








60 


216 


150 


5000 








75 


144 


150 


4500 








85 


144 


150 


4500 


B 8 


30 


0.5 


5 


352 


150 


10000 








15 


344 


150 


9500 








> 30 




Too broad \C) ll 




C 8 


10 


0.5 


5 


296 


150 


10000 








15 


304 


150 


9500 








30 


280 


150 


8500 








45 


232 


150 


6500 








60 


216 


150 


5000 








75 


176 


150 


4500 








85 


160 


150 


4000 


D 8 


1 


0.5 


5 


224 


80 


5000 








15 


216 


80 


5000 








30 


216 


72 


5000 








45 


200 


80 


5000 








60 


192 


80 


5000 








75 


160 


88 


5000 








85 


128 


88 


5000 


E 2 


10 


0.5 


5 


288 


150 


8500 








15 


296 


150 


8000 








30 


264 


150 


7000 








> 45 




Too broad [0 i] 




F 1 


10 


0.5 






Too broad I 




G 1 


1 


0.5 




144 


96 


3500 



Table 3 

Half Line Widths 



"The ratio of the kinetic energy to E initially deposited. 

''Angle between the line of sight and the jet direction. 

'^Half line widths of the [Fe n]-blend (near 5200A) for SN 1998bw. The observed value is 380 A. 
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Fig. 1. — The p — T conditions of individual test particles at their temperature maximum (where Tg = T(K)/10^). Fpr 
model C, the open circles denote those along the z-axis, and the filled circles denote those along the r-axis. The open 
squares denote those of model F (spherical model). 
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Fig. 2. — The isotopic composition of the ejecta of model C in the direction of the jet (upper panel) and perpendicular to 

the jet (lower panel). The ordinate indicates the initial spherical Lagrangian coordinate (M^) of the test particles (lower 
scale), and the final expansion velocities {V) of those particles (upper scale). 
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Vr [cm/s] 

Fig. 3. — The 2D distribution of '""^Ni (open circles) and "^^O (dots) of model C in the homologous expansion phase. 
Open circles and dots denote test particles in which the mass fraction of ^^Ni and ^^O, respectively, exceeds 0.1. The 
lines are density contours at the level of 0.5 (solid), 0.3 (dashed), 0.1 (dash-dotted), and 0.01 (dotted) of the max density, 
respectively. 
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Fig. 4. The profiles of the [Fe ii] feature (upper panels) and of [O i] 6300, 6363 A (lower panels) for model C viewed at 
15° from the jet direction (left panels; thick lines) and for model F (right panels). The observed lines at a SN rest-frame 
epoch of 216 days are also plotted for comparison (thin lines, Patat et al. 2001). The intensities of the strongest lines, 
normalized to O I] 6300.3A are: [Fe il] 5158.8A: 0.122; [Fe ll] 5220.lA: 0.026; [Fe ll] 5261.6A: 0.083; [Fe ll] 5273.3A: 0.039; 
[Fe II] 5333.6A: 0.060; [Fe III] 5270.4A: 0.032; [O l] 5577.3A: 0.022; and [O l] 6363.8A: 0.330. 
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Fig. 6. — The profiles of the [Fe ii] feature for hypernove models viewed at 5°. Each panel shows model A (upper-left), 
B (upper-right), C (lower-left), and E (lower-right), respectively. 
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Fig. 7. — The profiles of the [Fe ii] feature (upper) and the [O i] (lower) for normal energetic models. Those of model D 
viewed at 15° (left), 75° (center), and model G are shown. 



